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Abstract 

CN I During de Sitter inflation massless particles of minimally coupled scalar fields acquire a mass 

and a decay width thereby becoming quasiparticles. For bare massless particles non-perturbative 

CS| ■ infrared radiative corrections lead to a self-consistent generation of mass, for a quartic self in- 

Ph . teraction M oc X^H, and for a cubic self-interaction the mass is induced by the formation of a 

. ^ I non-perturbative condensate leading to M oc XsHs . These radiatively generated masses restore 

de Sitter invariance and result in anomalous scaling dimensions of superhorizon fluctuations. We 
introduce a generalization of the non-perturbative Wigner-Weisskopf method to obtain the time 
evolution of quantum states that include the self-consistent generation of mass and regulate the 

i-C . infrared behavior. The infrared divergences are manifest as poles in A = M^/3-ff^ in the single 

I I particle self-energies, leading to a re-arrangement of the perturbative series non- analytic in the 

rr^l couplings. A set of simple rules that yield the leading order infrared contributions to the decay 

,S^ ' width are obtained and implemented. The lack of kinematic thresholds entail that all particle 

states acquire a decay width, dominated by the emission and absorption of superhorizon quanta 

CN . CSC (X/H)^'^ [H/kph{r))]^ ; A [H/kph{r))]^ for cubic and quartic couplings respectively to leading or- 

^ys i der in M/H. The decay of single particle quantum states hastens as their wavevectors cross the 

^P I Hubble radius and their width is related to the highly squeezed limit of the bi- or tri-spectrum of 

scalar fluctuations respectively. 
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I. INTRODUCTION 



Quantum fluctuations during inflation seed the inhomogeneities which are manifest as 
anisotropies in the cosmic microwave background and are responsible for large scale struc- 
ture formation. In its simplest inception the inflationary stage can be effectively described 
as a quasi-deSitter space time. Early studies |l|-l6| revealed that de Sitter space time features 
infrared instabilities and profuse particle production in interacting field theories. Infrared 
divergences in loop corrections to correlation functions hinder the reliability of the pertur- 
bative expansion |7H9J, led to the suggestion of an infrared instability of the vacuum |1Ch14|. 
and affect correlation functions during inflation |7|, |8|, Il5l-ll8| requiring a non-perturbative 
treatment. 

Cosmological expansion modifies the energy-uncertainty relation allowing "virtual" exci- 
tations to persist longer, leading to remarkable phenomena, which is stronger in de Sitter 



space time as clarified in ref . [19| . Particle production in a de Sitter background has been ar- 
gued to provide a dy namical "screening" mechanism that leads to relaxation of the cosmolog- 
ical constant |2Ch22| through back reaction, much like the production of particle-antiparticle 
pairs in a constant electric field. The possibility that back reaction from the production of 
virtual excitations may yield a dynamical mechanism of evolution of dark energy rekindled 
the interest on infrared effects in de Sitter space time. A body of work established that 
infrared and secular divergences are manifest in super-Hubble fluctuations during de Sitter 
(or nearly de Sitter) inflation 23|-|25| and also invalidate the semiclassical approximation 26|, 
thus a consistent program that provides a resummation of the perturbative expansion is 
required. One possible approach is furnished by the dynamical renormalization group 27 



which provides a non-perturbative resummation of the secular divergences and has been im- 
plemented in several studies in de Sitter space time|28| and suggests a dynamical generation 
of mass[26|. The generation of a mass through the build up of infrared fluctuations was 
originally anticipated in the seminal work of ref. 29| , and explored and extended in ref. 30(] , 
and more recently a self-consistent mechanism of mass generation for scalar fields through 



infrared fluctuations has been suggested [23|, l26|, I3ll434 



Another particular aspect of the rapid cosmological expansion is the lack of a global time- 
like killing vector which leads to remarkable physical effects in de Sitter space time, as it 
implies the lack of particle thresholds (a direct consequence of energy-momentum conserva- 
tion) and the decay of fields even in their own quanta 28|, |35| with the concomitant particle 



production, a result that was confirmed in ref. [12|, |36| and more recently investigated in 
ref. |37l . |38| for the case of heavy fields. For light scalar fields in de Sitter space time with 
mass M <ti H, it was shown in ref. |28| that the infrared enhancement of self-energy correc- 
tions is manifest as poles in A = M'^/SH'^ and that the most infrared singular contributions 
to the self-energy can be isolated systematically in an expansion in A akin to the e expansion 
in critical phenomena. A similar expansion was noticed in refs. 26|, |3l|, |38|, |39| . 



Most of the efforts towards understanding infrared effects in de Sitter (or quasi de Sitter) 
cosmology focus on correlation functions, and only recently the issue of the time evolution of 
the quantum states has began to be addressed. In ref. [40| the Wigner-Weisskopf method |4lL 
42| ubiquitous in quantum optics^] has been adapted and extended as a non-perturbative 
quantum field theory method in inflationary cosmology which allows to study the time 
evolution of quantum states. This method reveals how quantum states decay in time and 
it has been shown to be equivalent to the dynamical renormalization group in Minkowski 
space time [40]. 



Motivation and results: 

There are at least two reasons to study the evolution of quantum states by implementing a 
non-perturbative resummation method: i) questions of instability of vacuum or excited states 
cannot be directly addressed by studying correlation functions perturbatively, ii) infrared 
and secular divergences in correlation functions require a resummation of the perturbative 



series. Whereas the dynamical renormalization group |26l l27l . |38| provides a resummation 
scheme in some cases, alternative non-perturbative resummation methods may offer novel 
insights and undoubtedly will be a welcome addition to the non-perturbative tools to study 
dynamical phenomena in cosmology. 

The main observation is that in the interaction picture field operators feature the free 
field time dependence and all the interaction effects are contained in the time evolution of 
states, therefore a method that provides a non-perturbative resummation scheme for the 
time evolution of states may provide an equivalent resummation framework for correlation 
functions by saturating the intermediate states with the time evolved states obtained from 
the non-perturbative time evolution. 

In this article we combine the expansion in A advocated in ref.[28|] with the Wigner- 



Weisskopf method introduced in ref . [40| to study the nature of the single particle excitations 
during de Sitter inflation. 

We focus our study on massless minimally coupled scalar field theories with typical cubic 
or quartic interactions and find that, similarly to finite temperature field theory, these ex- 
citations become quasiparticles acquiring a self- consistent mass that regulates the infrared 
as found in refs. 26|, |29|, |38| but also a decay width. 



In section II we study the generation of radiatively induced mass through the build up of 
infrared effects in a self-consistent manner. In the case of cubic self-interaction vertex the 
strong infrared behavior leads to the formation of a condensate that reveals that the theory 
is driven to a new stable minimum by radiative corrections, the expectation value of the field 
in this state yields a self-consistently induced mass. In the case of a quartic self-interaction 
the resummation of tadpole-type diagrams which are infrared divergent in the massless 



theory lead to the self-consistent generation of a mass, confirming the results of refs. [26 



29| . This self-consistent mechanism regulates the infrared behavior and induces, radiative 
anomalous scaling dimensions of superhorizon fluctuations. The infrared singularities for 
massless particles are replaced by poles in A = A'P/3H^ with M being the radiatively 
generated mass. 

In section III we combine the self-consistent approach with a generalization of the non- 
perturbative Wigner-Weisskopf method to obtain the long time evolution of single particle 
quantum states. 

We find that single particle states decay via the emission and absorption of superhorizon 
quanta and obtain their decay "widths" both for super and sub horizon modes to leading 
order in the expansion in A. The self-consistent mass generation regulates the infrared 
behavior which is now manifest as poles in A and leads to a rearrangement of the perturbative 
expansion which is non-analytic in the couplings. The decay of quantum states hastens as 
their wavevectors cross the Hubble radius. We argue that the order of the poles in A reflect 
the number of superhorizon quanta emitted in the decay process and obtain a set of simple 
rules to extract the leading order contributions in A to the decay "widths". We provide 
an intepretation of the decay "width" of superhorizon modes in terms of a relation between 
the single particle self energy and the bispectrum (for cubic coupling) or tri-spectrum (for 
quartic coupling) of scalar fluctuations in a highly squeezed limit. The order of the pole in 



A reflects the number of squeezed sides in tlie bi-or tri-spectrum configuration respectively. 
Conclusions, comments and further questions are presented in section IV. An appendix 
is devoted to the calculation of the self-energy for a cubic coupling to leading and next to 
leading order in A. 



II. 



CONDENSATE AND SELF-CONSISTENT MASS GENERATION: 



We consider scalar field theories in a spatially fiat Friedmann-Robertson- Walker (FRW) 
cosmological spacetime with scale factor a{t). In comoving coordinates, the action is given 
by 



S 



d^xdta^{t){ -0^- 



^4K-^-K-4 



p = 3,4 



with 



7^ = 6( - + - 

a a^ 



(2.1) 



(2.2) 



being the Ricci scalar, ^ = 0, 1/6 correspond to minimal coupling and conformal coupling 
respectively. 

Specializing now to the de Sitter case with a{t) = e^*, it is convenient to pass to conformal 
time 1] = —e^^^/H with dr] = dt/a{t) and introduce a conformal rescaling of the fields 



a{t)(p{x,t) = x{x,v)- 



(2.3) 



The action becomes (after discarding surface terms that will not change the equations of 
motion) 



S 



1 



d'xdvi;- x" - (Vx)' - M\v) x' - ^[c{v)f'''\'' 



with primes denoting derivatives with respect to conformal time rj and 

C'iv) 



where for de Sitter spacetime 

Civ) = a{t{v)) = 



Civ) ' 



Hv 

In this case the effective time dependent mass is given by 



M\v) 



1 



M^ / 1 



(2.4) 



(2.5) 



(2.6) 



(2.7) 



in what follows we consider the case of minimal coupling to gravity, namely ^ = 0. The 
Heisenberg equations of motion for the spatial Fourier modes of wavevector k of the fields 
in the non-interacting (A = 0) theory are given by 



x'Kv) + 



k'-^h' 



i(^'-^)]x.K^) = o, 



(2i 



where 

"^ = I - 171 (2.9) 

We will choose Bunch-Davies vacuum conditions for which the two linearly independent 
solutions are given by 

gu{k;v) = ^t'^-^V^^Hi'\-kr]) (2.10) 

U{k;v) = lt-''-'-V^^Hi'\-kv)=g:{k;v) , (2.11) 

where H^ ' (z) are Hankel functions. Expanding the field operator in this basis in a comov- 
ing volume V 

X{x, ^) = -^ E ["kaAk; r^) e*-^ + a\ gHk; r^) e-'^'^] . (2.12) 

k 
The Bunch-Davies vacuum is defined so that 



a,HO) = 0, (2.13) 

In the Schroedinger picture the quantum states |^E'(^)) obey 



and the Fock states are obtained by applying creation operators al to the vacuum 



i^\m{r^)) = H{r^)\m{r^)) (2.14) 

where in an expanding cosmology the Hamiltonian H{ri) is generally a function of rj. Intro- 
ducing the time evolution operator U{r],r]Q) obeying 

'-^Uiv, Vo) = H{r]) U{v, Vo), U{vo, Vo) = 1, (2.15) 

the solution of the Schroedinger equation is |^(^)) = U{ri,rio) [^(r^o)). Writing the Hamilto- 
nian as H{ri) = HqIt]) + Hi{ri) with HqIt]) the non-interacting Hamiltonian, and introducing 
the time evolution operator of the free theory f/o(^)''?o) satisfying 

(2.16) 
the interaction picture states are defined as 

l*(^))/ = Ujiv,Vo)\'^iVo))i = U,-\v,Vomiv))- (2.17) 

where Ui{ri,rio) is the time evolution operator in the interaction picture obeying 

^Uj{v,Vo) = -tHi{r])Ui{v,Vo), Ui{vo,Vo) = 1 (2.18) 

ari 

where the interaction Hamiltonian in the interaction picture 

Hiiv) = U,\r],r]o)H,ir])Uoiv,Vo) (2.19) 



is given by 



Hiiv) 



A 



[x{x,v)Y ^^^ 



and X is the free field Heisenberg field operator in eq.f l2.12p . 
In perturbation theory 



(2.20) 



Ui{v,m) 



l-i 



dr]'Hi{r]') + 



no 



f2.211 



In the interaction picture operators evolve in time with the free Hamiltonian Hq whereas 
states evolve as in eqn. fl2.17p . 

A. The tadpole and the A expansion: 



The tadpole will play an important role in the mechanism of self-consistent mass gener- 
ation, it is given by 



{0\x\x,r])\0) 



9u{k,V)\ 



(27r 



in 7]^ 



dz 



\H^\z)t 



for the massless, minimally coupled case v = 3/2 and 

,2 2 



\Hi'\z)\ 



l + z' 



n 



(2.22) 



(2.23) 



in which case the integral features both the usual quadratic and logarithmic ultraviolet 
divergence, but also a logarithmic infrared divergence. For minimally coupled "light" fields 
with W^/H"^ < 1 it follows that 



2 



A 



MP 



+ 



and 



W\z) 



2 z^O 



2" T(u) 



vr 



.2 A 



(2.24) 



(2.25) 



thus A > regulates the infrared behavior of the tadpole. To isolate the infrared we intro- 
duce ultraviolet (Ap) and infrared {^p) cutoffs in physical momenta and write the integral 
(I222D as 



H dz 



\Hl'\z) 



H dz 



z' \H^P(z)\ + 



H dz 



tfp z 

H 



MHz) 



(2.26) 



yUp — )■ acts here as infrared cutoff for the first integral. The second integral is ultraviolet 
and infrared finite for finite /ip, Ap. We can set u = 3/2 in this integral and use (I2.23p . 
In the first integral we obtain the leading order contribution, namely the pole and leading 
logarithm, by using the small argument limit of the Hankel functions 02.25P and we find that 
eq. 02.26p yields after calculation. 



H dz o I rn/ X |2 2 

-z' Hi'\z) =- 
z ' vr 



+ ;^ + 7-2 + ln% + 0(A) 



2 A 2H^ 



H 



(2.27) 



where we have displayed the pole in A and the leading infrared logarithm. In the second 
integral in f l2.26p we set v = 3/2 and combining its result with (12.271 ) we find that the 
dependence on the infrared cutoff fip cancels in the limit /ip — )■ leading to the following 
final result for the tadpole 



(0|x'(x,r/)|0) 



^^2^2 






(2.28) 



where 7 is the Euler-Mascheroni constant. While the quadratic and logarithmic ultraviolet 
divergences are regularization scheme dependent, the pole in A arises from the infrared 
behavior and is independent of the regularization scheme. In p articular this pole coincides 
with that found in the expression for < 0^(x, t) > in refs. |26l. |28| . l3ll |39| . The ultraviolet 
divergences, in whichever renormalization scheme, require that the effective field theory 
be defined to contain renormalization counterterms in the bare effective Lagrangian, for 
the tadpole this counterterm is of the form xiv) J{v) S'^d J{n) is required to cancel the 
ultraviolet divergences. Thus, the renormalized tadpole 



l{r]) = {0\x\x,r])\0), 



1 



87r2r/2 A 
where the dots stand for higher order terms in A ^ 1. 



[1 + ...] 



(2.29) 



B. Self-consistent mass generation: 

1. \x^ theory: condensate formation 



In this theory radiative corrections induce an expectation value of the field in the 
"dressed" vacuum state, which up to first order in perturbation theory is given by 



Hv)) 



l + z- / -^ d^x x'(^,V) 



|o) 



and the expectation value of x to leading order in A is given by 

A n d7]' 



{0x{y,vm = 3t 



H j„„ ri' 



vo 



d'^x 



X{y,v),x{x,r]') {0\x {x,r]')\0) 



this expression is depicted in fig. ([I]). We find 

3 A r dri' 



{o\xiy,vm 

where 



IQti H u ./,,(, ^'^ 



P+ (..'\P- 



{riY+ (r/ 



'\P+ r^^/3- 



{ri'f- iv) 



dz 



z' \H^;\z) 



'^ 2 



(2.30) 



(2.31) 



(2.32) 



(2.33) 



In order to understand the nature of the infrared divergences, let us first consider the 
massless case, namely A = 0, for which we find 



(o|x(f,^)|o) 



3A 



IGn H u 7] 



In 



^0 



dz 

z 



z'\H^Hz)\\ 



(2.34) 




FIG. 1: Tadpole contribution to the expectation value (0|x(y,f?)|0). 



There are two sources of infrared singular physics in this expression, i) for u = 3/2 the 
z integral is infrared divergent (it is also ultraviolet divergent but this divergence can be 
canceled by the counterterm discussed above), ii) in the long time limit rj ^ 0" the loga- 
rithmic term is secular and entails that the (unsealed) expectation value grows in time: the 
factor —l/Hr] = a{t) reflects that the expectation value of the unsealed field (p = x/o-if) 
(12 .Sp would be constant were it not for the logarithmic term. The growth of the expectation 
value implies the formation of a condensate. For the massless case (0|x(y, ^)|0) — ln^[?7/?7o] 
if the integral is regulated with an infrared cutoff constant in comoving coordinates [26|. 

For M^ 7^ both infrared divergences are regulated, using the result (I2.29P and A ^ 1 
we find for rj/rjQ — )• 0, 



(o|x(^,r/)|o) = x(^) 



A 



— > 



[1 + 0(A) + ---]. 



The un-scaled field (f) acquires a constant expectation value asymptotically for r^/r/o — > 

1 ,^, ,^ .,^, A 



(2.35) 



mmrim 



mxiy.rim^ 



[1 + 0(A) + ---]. 



Since the field is acquiring an expectation value we shift the field and define 
X(f,r/) = ^(f,r7) + x(r/) ; (0|^(f,r7)|0) = , 
introducing this shift in the interaction Hamiltonian (I2.20p for p = 3 we find 



H, 



where 



leading to 



d?x 



1 M2 



1 M2 . 

7^2 2 if 2 






vj/- 



M2 



7^2 2 if 2 



3A2 



-'^^^^) 



1 + 0(A) 



(2.36) 



(2.37) 



(2.38) 



(2.39) 



(2.40) 



47r2ii2A2 

and we neglected terms that are constant and linear in \E' (see the discussion below). This 
suggests a mechanism of self- consistent mass generation, indeed interpreting M as the mass 
of the field, with A given by (I2.24p eqn. (I2.40p becomes a self-consistent condition with the 
solution 



M = H VS 



- A 
.27rii. 



1/3 



(2.41) 




FIG. 2: The tadpole and expectation value contribution to the self-energy in Xx^ 



This mass term is identified with a self-energy contribution depicted in fig. (j2]). 

A systematic implementation of this self-consistent mechanism can be formulated as 
follows: 

• As in renormalized perturbation theory, a perturbative expansion around the vacuum 
state with the correct mass including radiative corrections is achieved by adding a mass 
term to the non-interacting Lagrangian density and subtracting it in the interacting 
part as a mass counterterm, namely 



C 

Co 



Co + Cj 
1 



x"-i^xr-^.{^'^)x^ 



M2 



A 



2 if 2 ^2 



X 



-Hv] 



X 



(2.42) 
(2.43) 



Expand the field in Fourier modes in a finite spatial comoving volume V quantizing 
the modes with the usual periodic boundary conditions in this volume 



X{x,r]) 



1 



Y^ X{k, v) e 



ik-x 



Xik,ri) = af:g^{k;ri) + al^gl{k;ri) 
where the mode functions are given by (12.101) with u given by (12. 9p . 



(2.44) 
(2.45) 



Separate the zero mode x(0, rj) and normal order the full Hamiltonian in the interaction 
picture of the free massive field with mass M. The cubic interaction becomes 



2^ 



£xx^ix,v) = -i=^T—{^\y2-Xik,v)xi-k,v)- +VI{r]) x(0, r/) 
yVHri l-'t^ j 



+ X^ ■X{k,vi)x{,q,ri)x{.-k-q,vi) ■.+x^{Q,vi)) (2.46) 

where we have renormalized the tadpole by canceling the ultraviolet divergences with 
appropriate counterterms yielding the finite contribution X(?7) given by (12.291) . In the 
V^ — ;■ oo limit the expectation value of the field 



{x{x,r])) = X{V) 



1 



(X(0,r/)) 



(2.47) 



is now determined by the linear term \/VI{ri) x{0,v) ^^id we find from the result 

(X(0,.))^-^^^ (2.48) 

The presence of a vacuum expectation value of x depicted in fig. ([T]) indicates that 
an effective action for a long- wavelength component of the field features a linear term 
in X along with the induced mass term from this expectation value indicated in fig. 
([2]). Such a linear term indicates a new, radiatively induced minimum of the effective 
potential and requires that the field be shifted by the new vacuum expectation value, 
so that the one-point function (expectation value) of the shifted field vanishes in the 
new (correct) vacuum. Therefore, as usual, quantization around the condensed state 
is achieved by shifting the field by its expectation value, defining "^ as in eqn. (I2.37p . 
carrying out the Fourier transform and introducing the linear counterterm J^{ri) that 
ensures the vanishing of the expectation value of \E' one finds^ 



£--1 — — t: 



M^ 3A^ 



V1/2+ ^ : vl/3 : + {j{v) + jf- [l{v) + x\v)] }^ (2.49) 



Hr] I Hi] 



with : \E'^ := \1''^ — 3X{r]) where X{f]) is given by (I2.29P and we have neglected constant 
terms. The linear counterterm J7'(?7) is fixed by requiring that (0|\E'|0) = in the 
correct vacuum state systematically in perturbation theory. To one-loop order 

3A r 1 

J{v) + -Fr- Av) + x\v) =0. (2.50) 

H rj L J 

This choice of counterterm simultaneously ensures that (0|^|0) = and the cancel- 
lation of the matrix element of the interaction between the vacuum and the single 
particle state with zero momentum (l(j|if/|0) = 0, thus ensuring harmonic perturba- 
tions around the new vacuum state. Cancellation of the quadratic term in Ci yields 
the self-consistent gap equation leading to the radiatively generated mass given by eqn. 
(12.411) . however, we leave the mass counterterm in Cj anticipating the possibility of 
a further one loop contribution from the non-local self energy diagram studied in the 
next section. 

It is important to highlight that renormalizing the short distance divergence with 
different renormalization schemes in < ^^ > does not affect the leading order pole in 
A, a purely infrared effect, as different subtraction schemes differ by finite constants. 
Thus the self-generated mass term is a genuine radiative infrared effect. 

The fact that one (zero momentum) mode of the x fi^ld becomes macro scopically occu- 
pied, with {x^iV))"^ oc V^ in the V^ — )■ oo limit (see eqn. (I2.48P ) signals the emergence of a 
condensate as a consequence of radiative corrections. As mentioned above the rj dependence 
of this condensate entails that the zero momentum Fourier mode of the unsealed field 
acquires a constant and macro scopically large expectation value. In Minkowski space time 
as well as in de Sitter space time with conformally coupled fields, the ultraviolet subtracted 



^ A linear term j^l i X^ is cancelled between the free and counterterm parts of the Lagrangian density. 

10 



one-point function vanishes in the niassless limit, thus both the self-generated mass term 
and the condensate are genuine radiative infrared effects in de Sitter space time and fields 
minimally coupled to gravity. 

The self-consistent mass fl2.4ip corresponds to 



A 



271 H. 



hence consistency with the approximation A ^ 1 requires that X/H ^ 1. 
Furthermore, it follows from the above analysis that 



(02(x,r7)) 



1 



C'iv) 



r{v)+Av) 



ocH' 



H 

T 



2/3 



namely a de Sitter invariant result but non-analytic in the coupling. 



(2.51) 



(2.52) 



2. Ax^ theory: 

For p = A the interaction Lagrangian density is £/ = — Ax^ the mechanism of self- 
consistent generation of mass via a tadpole contribution is implemented as follows: add a 
mass term to Cq and subtract it as a counterterm in Cj leading to 

X'-Xx'. (2.53) 



Ci 



2 if 2 ^2 

Requiring that the counterterm cancels the contribution from the tadpole depicted in fig. 
dH]) to the two-point function yields the condition 

M2 



12X{x\x,r])) 



(2.54) 



which upon renormalization by subtracting the ultraviolet divergences and using the result 
(I2.29P leads to the self-consistent condition 

M2 12 a 



with the solution 



M 



ff2 ^2 3^2 ^2 ^ 

9 



L27r2J 



X^H. 



(2.55) 



(2.56) 



The coupling dependence of this result is in agreement with those of refs. |2| 
The self-consistent mass ( 12.56^ leads to 



A 



A 

2^2 J 



321,13^ 

(2.57) 



and again, consistency with the approximation A ^ 1 requires that A ^ 1. 

After the usual ultraviolet renormalization the equal time two point function of the un- 
sealed scalar field (see eqn. (12. 3p ) is de Sitter invariant, non- analytic in the coupling and 
given by 

(</)^(r^)) oc ^ (2.58) 

which agrees with the result in refs. 29|, |4^ . 



11 




FIG. 3: Tadpole contribution to the self-energy in Ax 



C. Normal ordering vs. self-consistent mass generation: 

A simple approach would be to normal order the interaction Hamiltonian, which amounts 
to subtracting tadpole contributions thereby avoiding the problem of infrared and ultraviolet 
divergences. However, normal ordering requires specifying a vacuum state, as can be seen 
from the identity 

: Xkiv)X-k{v) ■■= Xk{v)X-kiv) - loAk; V)\^ (2-59) 

in particular normal ordering in the bare massless theory corresponds to z/ = 3/2 leading to 
the infrared divergence of the tadpole, whereas normal ordering in the theory with a mass 
term in the non-interacting Lagrangian implies u ~ 3/2 — A and the infrared divergence 
is regularized, for conformally coupled massless particles z/ = 1/2, and the tadpole only 
features the ultraviolet divergence of Minkowski space time but no infrared divergences. 

In the case of the cubic vertex for minimally coupled massless fields naive normal ordering 
neglects the fact that there is an infrared and secularly divergent radiatively induced vacuum 
expectation value, thus the vacuum state evolves to a state in which the one point function 
is non-vanishing. Furthermore, this evolution of the expectation value is leading to the 
generation of a radiatively induced mass, which in turn regulates the infrared behavior. 
Similarly for the quartic vertex, the contribution from the tadpole and the self-consistent 
resummation entails that the vacuum state evolves on to another state in which the infrared 
divergences are regulated by the self-consistent mass. 

Furthermore and important consequence is that the original massless theory breaks 
the underlying de Sitter invariance because of the infrared effects [4, |5|, whereas the self- 
consistent generation of mass restores de Sitter invariance by manifestly regularizing the in- 
frared and leading to a time independent equal time expectation value (0^(^)) = constant, 
as befits de Sitter invariance. Thus the mechanism of self-consistent mass generation re- 
markably restores de Sitter invariance in the case of bare minimally coupled massless fields. 

Last but not least, whereas normal ordering discards the tadpole divergences, the infrared 
divergences emerge in non-local self-energy contributions at one loop in the case of the cubic 
vertex or two loops in the case of the quartic vertex, this will be seen in detail below. 

Through the self-consistent mass generation, the infrared divergences become poles in 
A and relieve the large logarithms at the expense of a rearrangement of the perturbative 
expansion non-analytic in the coupling. 

Thus normal ordering from the outset in the massless theory misses the important physics 
associated with the buildup of infrared effects, the fact that time evolution is rearranging the 
vacuum (and many particle states) in a manner that ultimately self-consistently regulates the 
infrared behavior. This situation is similar to the results from the stochastic approach |29l. 



30. 
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D. Anomalous dimensions from self-consistent masses 



While the emergence of self-consistent masses has been previously recognized [26|, |29|, |31, 



32|, |38|, it is noteworthy that this infrared mechanism of mass generation leads to anomalous 
scaling dimension of the two point function in the superhorizon limit. In the bare mass- 
less theory the two point function scales as 1/k^ in this limit but after the self-consistent 
generation of mass these now behave as 

ixdv) X.dv)) ^ ^ ■ (2.60) 

where A is given by eqns. fl2.51ll2.57p for cubic and quartic couplings respectively. We 
emphasize that the anomalous scaling dimension A is a non-perturbative result from the 
build-up of infrared fluctuations, much in the same way as anomalous scaling dimensions 
in the theory of critical phenomena. Furthermore, in the case of the cubic coupling, it is a 
result of the formation of the condensate, which is radiatively induced and a consequence of 
the infrared divergences of the massless theory. The emergence of anomalous dimensions as 
a consequence of the infrared generated mass has also been recognized in ref . |25| . 

In the next section we show how the self-consistent method is systematically implemented 
within the Wigner-Weisskopf theory and also study how quantum states acquire a decay 
width. 



III. WIGNER-WEISSKOPF THEORY : 

A. Transition amplitudes and probability 

In anticipation of the non-perturbative Wigner-Weisskopf theory in de Sitter cosmology 
and to identify corrections to masses and widths of the states, let us consider the example 
of two interacting scalar fields, x? ^ after conformal rescaling, with interaction Hamiltonian 

Hi = ^- [d^x x{x,v) S\x,7]). (3.1) 

-Ht] J 

The case (5 = x is obtained straightforwardly. Using the expansion of the scalar field x given 
by fl2.12p and expanding the field 6 as 

Si^^v) = ^J2 [h ^^(k; v) e^'-' + h\ v^{k; v) e^^^-^j . (3.2) 



where 



with 



uAk,v) = t; i''^^^V^^H^^\-kri) (3.3) 



.^^^§. (3.4) 



The transition amplitude x ~^ 25 is given by 



~^x^-25(kj]ri) = -—j=l ^g^{k;rii)u^{p;rii)u^{q;rii) ; q=\p + k\ (3.5) 
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and the total transition probability is 

rv p 

T'^^2s{k; T]) = dr]2 dr]i T,{k ■,T]i,r]2) (3.6) 

•Jrin J rin 



'170 Jrjo 



wliere^ 

^{k;7]i,r]2) = j^ / j^Uj,{p,7]i)uTj{q,r]i)uTj{p,r]2)UTj{q,7]2) ; q = \p+k\ 

(3.7) 
with the property that 

^ik;ri2,rii) = T.*{k;r]i,r]2). (3.8) 

Introducing the identity 1 = Q{r]2 — i]i) + Q{r]i — ri2) in the (conformal) time integrals and 
using fl3.8p we find 

rv rv2 p -, 

V^^2s{k;r]) =2 dr]2 dr]i Re S(A; ; 771,772) (3.9) 

from which we identify the transition rate 

r(77) = ^V^^25{k;r]) = 2 P dr]' Re[E{k;r],r]')] (3.10) 

The result can be extrapolated to the self-interacting Ax^ case simply by replacing the 
mode functions u-uik^rj) — > g^{k,r]) and including the corresponding combinatoric factor. 

In Minkowski space-time 77 — )■ t, if the kinematics of the transition is allowed, namely 
energy-momentum conservation holds in the process, the transition is to on-shell states and 
the transition probability features a secular growth linear in time at long time, in which 
case the transition rate becomes a constant. This is the result from Fermi's Golden rule. If, 
on the other hand energy-momentum conservation is not fulfilled, the probability becomes 
constant at asymptotically long times, with a vanishing transition rate, describing virtual 
processes that contribute to wave function renormalization. A true decay of the quantum 
state is therefore refiected in a secular growth of the transition probability and a transition 
rate that either remains constant or grows at asymptotically long time. In de Sitter space 
time the lack of a global time-like Killing vector implies the absence of kinematic thresholds 
and the lack of defined "on-shell" states. As discussed earlier in ref . |28l . |35| quanta of a single 
field can decay into other quanta of the same field, and more recently|40|] a generalization 
of the non-perturbative Wigner-Weisskopf method to cosmology was introduced to study 
explicitly the decay of quantum states. 

B. Wigner-Weisskopf theory in deSitter space time: 

In order to make the discussion self-contained, we highlight the main aspects of the 
Wigner-Weisskopf non-perturbative approach to study the decay of quantum states pertinent 
to the self-consistent description of mass generation discussed in the previous sections. For 



^ A factor 1/2! accounts for Bose symmetry of the two particle final state. 
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a more thorough discussion and comparison to results in Minkowski space time the reader 
is referred to ref.|40|. Expanding the interaction picture state \'^{t]))j in Fock states \n) 
obtained as usual by applying the creation operators on to the (bare) vacuum state (here 
taken to be the Bunch-Davies vacuum) as 

\^{v))i = Y.Cn{v)\n) (3.11) 

n 

the evolution of the state in the interaction picture given by eqn. 02.17P yields 

i^\^iv))i = Hjiv)\<ilir^))j (3.12) 

which in terms of the coefficients Cn{i]) become 

dCniv) 



di] 



J2Cn.iv){n\Hjiv)\m) , (3.13) 



m 



it is convenient to separate the diagonal matrix elements, that represent local contributions 
from those that represent transitions and are associated with non-local self-energy correc- 



tions^, writing 

^^M = -^C^(r^){n\HJir^)\n) - i J^ CUv){n\Hj{v)\m) . (3.14) 

Although this equation is exact, it yields an infinite hierarchy of simultaneous equations 
when the Hilbert space of states \n) is infinite dimensional. However, progress is made by 
considering the transition between states connected by the interaction Hamiltonian at a 
given order in Hj: consider the case when one state, say \A) couples to a set of states \k), 
which couple back to \A) via Hj, to lowest order in the interaction the system of equation 
closes in the form 

^^ = -^(A|i^,(^)l^)C^(^)-^E(^l^^(^)l'^)C'«(^) (3-15) 

^^ = -tCAiv){n\HM\A) (3.16) 

where the J2k=/=a ^^ '^^^^ ^^^ ^^^ intermediate states coupled to \A) via Hj representing 
transitions. 

Consider the initial value problem in which at time t] = r]Q the state of the system is 
given by |\l'(r7 = ?7o)) = \A) so that 

CAivo) = l ; C,Mr] = r]o) = 0, (3.17) 

solving (13.161) and introducing the solution into (13.151) we find 

CM = -^ [\K\Hj{r]')\A)CAW)dv' (3.18) 

^^M = -^{A\HJ{rJ)\A) CAv) - [' ^a{v,v') CaW) dv' (3.19) 

dv Jvo 



In ref. 40[ the diagonal matrix elements, hence the local contributions were not included. 
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where^ 

^a{v,v') = J2{A\Hjir^)\K){K\Hj{r^')\A) . (3.20) 

Equation (3.19) makes manifest that J^Aii],"!]') is the retarded self-energy, since rj > r( . In 
ref. 40(] the correspondence between the Wigner-Weisskopf and the Dyson resummation of 



the propagator in Minkowski space-time is explained in detail and shown that the self-energy 
that enters in the evolution equation for the amplitudes is the retarded one, resulting in that 
the propagator features a pole in the upper half of the complex frequency plane (see eqns. 
(2.12), (2.30) in ref.|40|). The retarded self-energy is the one that enters in the description 
of an initial value problem as befits the evolution of the amplitudes from an initial state. 

In eqn. 03.16P we have not included the diagonal term as in 03.15p ^. it is clear from 03.18P 
that with the initial condition (I3.17P the amplitude of C^ is of 0{Hi) therefore a diagonal 
term would effectively lead to higher order contributions to (13.190 . The integro-differential 
equation (I3.19P with memory yields a non-perturbative solution for the time evolution of the 
amplitudes and probabilities, which simplifies in the case of weak couplings. In perturbation 
theory the time evolution of Ca{v) determined by eqn. (I3.19P is slow in the sense that the 
time scale is determined by a weak coupling kernel S^, hence an approximation in terms of 
an expansion in derivatives of Ca emerges as follows: introduce 

W{v,v')= r ^Aiv,v")df (3.21) 

so that 

^A{v,v') = ^,W{r],r]'), W{r],r]o) = 0. (3.22) 

Integrating by parts in eq. (l3.19p we obtain 

' S^(r/, V) CAiv') di = Wiv, T]) CAiv) - r Wiv, V') -^CaW) dv'- (3.23) 

The second term on the right hand side is formally of higher order in Hj, integrating by 



parts successively yields a systematic approximation scheme as discussed in ref. |40 
Therefore to leading order in the interaction we find 

Ca{v) = e-^^l^^"^''^'^'"', W{v',v') = ^{A\HJ{r^')\A) + T S^(V, V')rfV' . (3.24) 

This expression has a clear and simple intepretation in Minkowski space time by replacing 
r],r]' — !■ t,t' with the states \n) being eigenstates of the unperturbed Hamiltonian Hq, and 
Hj{t) = e'^°^Hi{Q)e-'^°^: the contribution {A\Hi{t')\A) = {A\Hj{0)\A) = 5E^)^^ is simply 

^ In ref.[40| it is proven that in Minkowski space-time the self-energy in the single particle propagator is 

given by iS. 
^ These diagonal terms represent local self-energy insertions in the propagators of the intermediate states, 

hence higher orders in the perturbative expansion. 
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the first order shift in the energy in elementary quantum mechanics and taking the long 
time limit t — > oo (with a convergence factor ie) 



i— >oo 



^Ait,t')dt' 



^m!ML^uEf + ^ 



KJ^A 



E^-Ea + ie 



(3.25) 



thus the imaginary part of the time integral yields the second order energy shift and the real 
part yields half of the decay rate a la Fermi's golden rule. Inserting this result in (I3.24p 
yields, in Minkowski space-time 



CA{t) 



^-iSEAt g-ri/2 



6E, 



6E^j^ + 5Ef 



(3.26) 



A more careful treatment of the long time limit also exhibits the wave function 
renormalization[40| and establishes the equivalence with the dynamical renormalization 
group 271] in Minkowski space time. 

Motivated by this interpretation we introduce the real quantities Sa{i]) ] ^a{v) ^s 






in terms of which 



c^(^) = e-^/;^^('7'w,-l/;r.(.'w 



(3.27) 



(3.28) 



When the state A is a single particle state, radiative corrections to the mass are extracted 
from Ea and 

r^(r^) = -|^ln 



\CA{ri)\^ 



(3.29) 



is identified as a (conformal) time dependent decay rate. We see from (13.291) that '^{jf) is 
exactly the same as expression (I3.10p . 

In Minkowski space-time Ea corresponds to the self-energy correction to the mass of 
the particle J40|, |43| and the program of renormalized perturbation theory begins by writing 
the free field part of the Lagrangian in terms of the renormalized mass and introducing a 
counterterm in the interaction Lagrangian so that it cancels the radiative corrections to the 
mass from the self-energy. Namely the counterterm in the interaction Lagrangian is fixed 
by requiring that for the single particle state \A) = |1^), in the long time limit rj' — )■ 0~ 



^idv') = {h\Hi{v')\W + I im 



Si(fc;?7',?7") dr]" = 



(3.30) 



We will implement the same strategy to obtain the self-consistent radiatively generated mass 
in de Sitter space time where equation (I3.30p will determine the self- consistent condition for 
the mass. In Minkowski space time, the condition (I3.30p is tantamount to requiring that 



the (real part of the) pole in the propagator be at the physical mass [40 
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C. Cubic vertex: 

We now consider the cubic interaction given by the Lagrangian density fl2.43p . In the 
interaction picture the fields are expanded as in (12.121) . to carry out the perturbative expan- 
sion in the state in which the expectation value of the field vanishes, the field x is shifted by 
its vacuum expectation value as in f l2.37p . and include the tadpole counterterm given by the 
last term in eqn. f l2.49p . thus replacing the interaction Hamiltonian by that obtained from 
fl2.49p and treating the quadratic term ^^ as a mass counterterm, namely from eqn. fl2.49p 

HM= Id'Al^,"^" -^^^'^-f^(^) + |^[^(^) + X^(^)l 1^^ (3-31) 

The counterterm J{ri) is required to cancel the linear term in ^ in f l3.3ip . leading to the 
relation ( I2.50p . simultaneously ensuring that 

(0|^|0) = ; (l,^=ol^/(^)|0) = 0. (3.32) 

We focus on the time evolution of single particle states of the field Xi namely |lg). The 
mass counterterm contributes only to the diagonal matrix element 

{h\Hi{Tl)\h) = Jp^\9u{Kti)\ ■ -^ = 4^2ff2A2 - ff^' (3.33) 

The interaction Hamiltonian connects the state |1^) to an intermediate state with two 
particles |lp';l^_-,) and also to the state |lg; l^-; Ip-; 1^_-,), the first state is connected back 
to |1^) by Hj in the self-energy contribution depicted in fig.(jl]-a) whereas the second state 
contributes to the vacuum disconnected diagram displayed in the same figure (jH-b). The 
vacuum diagram is subtracted out consistently in perturbation theory by redefining the 
dressed states constructed out of the dressed vacuum. See discussion in ref. kc 



<^ ^ 



(a) (b) 

FIG. 4: Lowest order contributions to the evolution of the single particle states. Fig. (a): irre- 
ducible self energy, fig. (b): vacuum diagram. 



The irreducible self-energy diagram (a) is given by 

io\2 p j3„ ^ ^ 

S(A;,r/,r/') = ^ g*{k;r])g^{k;r]') / ——^ 9u{p;v) 9*u{p;v') 9u{\k - p\;v) 9l{\k ~ p\;v') 

(3.34) 
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with g^(k;r)) given by eqns. (12.101) . We recognize eqn. (13.340 as the kernel in the transi- 
tion probabihty (I3.6|3.7p with the only difference being the combinatoric factors and the 
replacement M17 — > g^. Furthermore with G{k,'r],'r]') = g*{k;ri)g^{k;r]') it follows that 

^{k,ri,ri') dri' oc I j^^^^ j -^^ G{k,ri,ri') G{p,ri,ri') G{\k - p\,ri,ri') (3.35) 



vo 



(2vr) 



Ht]' 



and the imaginary part of the r^'-integral is proportional to the bispectrum of the scalar 
field (45|. an observation that will become important in the interpretation of the final result 
below. 

To simplify notation we introduce 



and carry out the angular integral by changing variables to 

1 
P +p'^ + 2kp cos(6') ' ^ d cos{6) 

in terms of which the self-energy is given by 



kp 



dq 



(3.36) 



(3.37) 



J:{k;ri;ri' 



9A2 



9t{k;v)9uik;ri') I[k,ri,ri'] 



where 



I[k,r],r]'] 



32 H^k 

Hyip^rj-^if]') J^^{k + p]ri7]') - J'^dk - p\;r] r]') 



and[4e 

Fu{q]r]T]') 



V.y{q]T]]T]')qdq 



rW, 



r(2) 



pdp 



(2), 



(3.38) 
(3.39) 

r(i) 



-gV) Hl'>{-q7^) Hl'M-qv') ' {-QV) m'\-qi) K'l.^-qr,) 



rf — 1]'"^ 



(3.40) 

The integrals feature infrared divergences at p — t- and \k — p\ — )■ a result that follows 
from the identities 



n^{p;v;v') 

^u{q]r]]T]') 



p-^O 



(j-s>0 



rfz/)i2 / 4 



TC 



rji]' 



p 



-2u 



r(i/)r(i/-i) / 4 



2-2;^ 



2'k'^ \r]r]' 

The integrals near the region p ~ and p c:^ k yield simple poles in 



A = --Z/. 
2 



(3.41) 
(3.42) 

(3.43) 



In the appendix we carry out the integrals keeping the poles in A and the leading infrared 
logarithmic contribution. The final result for I[k,ri,ri'] in eqn. (I3.38P up to the leading and 
next to leading order in A is given by eqn. flA17p . 

The first term in the first line of eqn. (IA17P is the leading order contribution in the 
A — )■ limit and encodes the leading infrared behavior. The second term in the first line 
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is recognized as the contribution from a massless particle confornially coupled to gravity 
and yields the leading contribution for modes deep inside the Hubble radius, since for these 
modes the behavior of the Hankel functions is described by the simple Minkowski space time 
plane waves 40| . These two terms yield 



j:{k,r],r]') 



/\\2 



9y \g,ik,r])\'\g,ik,r]')\ 



2HH^ 



+ 



9A2 



8HH 



2^2 



[rjrj'Y A 
9l{k]ri)g^{k]r]') 



[fc^r^r]'] 



-ik{rj—ri') 



Hjir]') 



V\ 



1 



TT 



T] — rj' 



— Kv - v') 



. (3.44) 



Although the term in the second line is subleading in A, we have included it because it 
yields the leading contribution from subhorizon modes which are completely described by 
a conformally coupled massless field and contain the short distance ultraviolet divergence 
associated with the usual mass renormalization. 

From the result eqn. (IA20p . it follows that in the integral 



rv 

W{ri,ri) = / Tj{k,ri,ri') drj' 



(3.45) 



the second term in (I3.44p yields the contribution found in ref.[40|] in the case of massless 
conformally coupled particles 



Wcciv^ri) 



9\^\gu{k,r,)\ 



TC 



+ i In (e) 



(3.46) 



The imaginary part is an ultraviolet divergent contribution to the mass, which is canceled by 
the mass counterterm (I3.33P whereas the real part gives the contribution to the decay width 
found in ref. 40| for conformally coupled massless particles, which, however is subleading in 
A. The first term is purely real and does not contribute to the mass thus the results 03.33P 
and (I3.46P combined with the condition (I3.30p leads to 



jj2 



6M^ 



+ 



9X' 



L H^ gnm^ 



In (e) 







Absorbing the ultraviolet divergence in a renormalization of the mass M, 
condition becomes 



2 
ren 



(3.47) 



the self-consistent 



Mj 



3X' 



m 47r2 if 2^2 



which yields 



A. 



A 



V2txH\ 



Mr. 



H V3 



A 



27cH 



1/3 



(3.48) 



(3.49) 



confirming the result (I2.4ip in the previous section but now including the short distance 
renormalization. In what follows we will refer to A^en? Mren simply as A, M with the im- 
plicity understanding of renormalized quantities. 

Although the self-consistent mass (I3.49P is similar to the result in ref. [3l|, we note that our 
result reveals that the irreducible self-energy features a single pole in A, in agreement with 
the perturbative study in ref. 39(] and the "infrared counting" of ref. [26| but in disagreement 
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with the results of the plane wave ansatz proposed in ref. SjJ^ and that the self-consistent 
mass arises from the tadpole (expectation value), it is the result of the formation of the 
infrared condensate, and not from the irreducible self-energy which is real to leading order 
in A. These differences notwithstanding, the self-consistently generated mass features the 
coupling and H dependence in agreement with the results of ref. [3l| with a slightly different 
factor. 

Superhorizon modes: For —kr], —krj' <^ 1 the first line in eqn. f l3.44p dominates and 



S(fc,r/,r/') 



9\^e \H^u\-kri)\'' \Hl'\-kr]')\ 



/\|2 



32if2A [-kr]^-^ [-k7]'Y-'^ 

Therefore requiring the condition fl3.47p we find from the results 



(3.50) 



W{ri,ri) 



T{r]) 9\^k \m'{z)\ 



(1)^ 



32H^A 



vl-A 



'° \hI^Hz' 



r'U-A 



dz' 



(3.51) 



where z = —kri,ZQ = —krjQ. The integral in (|3.51D is remarkably similar to an expression that 
emerges in the two point correlator in refs. [ll|, 38] upon the analytic continuation u — )■ —in 
where /x = [^ ~ f] ^ > ^'^d real in those references. However for u = —ifi as in ll|, |38| the 
limits of the integral can be taken zq — ?► oo, z ^ 0, whereas the lower limit cannot be taken 
to vanish for z/ ~ 3/2 — A and real because of the strong infrared divergence, a consequence 
of light masses. This is an important difference with the case studied in these references 
that prevents a meaningful comparison. 

We can now obtain W{ri, rf) in eqn. (I3.24p . since S is real to leading order in A, requiring 
that the tadpole cancels the imaginary part of W and carrying out the remaining integrals 
in the limit 2; ^ 1 ; ^o — 1 (when the particular mode crosses the Hubble radius) we finally 
find 

i,(r/) ^ e-^^-^") . 7(-A:r/) = — ^^[-A:r/]~6(i-^). (3.52) 



Ci 



167r2iJ2/\,' 



Using the result (I3.49P we find 



l{-kri) = ^ 



2t^H 



H 



kphiv) 



6(1-A) 



+ 



(3.53) 



The dots stand for higher order terms suppressed with respect to the leading term by at 
least a factor [A/if] . Thus we see that the self-consistent mass generation leads to a 
rearrangement of the perturbative expansion non-analytic in the coupling A. 

The strong suppression of the single particle amplitude for subhorizon modes is a con- 
sequence of emission and absorption of soft superhorizon quanta. The dependence on the 
wavevector oc 1/k^ has a simple interpretation in terms of the relation of T,{k,rj,rj') with 
the bi-spectrum|45l] highlighted by eqn. (13.350 : the pole in A arises from the integration in 
a small band \p\ < fi -^ oi the highly squeezed triangle for the bispectrum configuration 
displayed in fig. ([5]). The extra power 6 A in (I3.53P reflects the anomalous dimension. 



^ The author has not been able to understand the source of the discrepancy, the plane wave ansatz proposed 
in ref. 



3l( does not seem to include the non-locality in conformal time explicit in the result (I3.44p and 



that of ref. 



391 
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\p\<^l^^ 



k — V 



FIG. 5: Triangle of momenta for the bispectrum (see eqn. (j3.35p ) integration over |p| < /i 
corresponds to the highly squeezed limit and yields the pole in A. 



From the result for the self-energy obtained in the appendix it is clear to next order in 
A oc [A/iJ] that there are logarithmic corrections oc A ln[— A;?]] to j^—kr]) thus these cor- 
rections will become relevant when —kr]^ e^l^ when the single particle amplitude is already 
strongly suppressed. These corrections may be included in the self-energy and the non- 
perturbative Wigner-Weisskopf resummation automatically exponentiates these potentially 
secular terms into a correction to the exponent '-^{—ki]) in ( 13.521) . Clearly extra logarithms 
in higher orders in A only modify the decay law but their "secularity" is automatically 
resummed into the overall decay function '-/{—kr]). 

Subhorizon modes: For —krjQ ^ —kf], —kt]' ^ 1 close inspection of the self-energy 
contributions given by eqn. flA17|) shows that only the first two terms displayed in the first 
line of the equation are dominant, however the contribution from V[k,7],7]'] is suppressed 
with respect the term 1/A -|- ln[— /c?]] in the bracket. The second term in the first line is 



given by eqn. flA20p . Furthermore, in the subhorizon regime 



9t{k;r])g^{k;r]') 



2k 



(3.54) 



after absorbing the short distance divergence in a renormalization of the mass and imposing 
the condition (I3.48P we find 



W{7],7]) 



9\H 



y3-A 



+ 



ttA 
4^2 



; z = -krj = kph{r])/H . 



(3.55) 



The final integral is performed in the limit —krjQ ^ —kr] ^ 1 and using (I3.49P leads to the 
final result for sub-horizon modes 



Cikiv) 



-yiv) 



liv) = 7 



9 / A 



4 \2nH 



H 



kphiv) 



2-A 



TT / A Y kphiv) 
2 \2nHJ H 



+ •■■ (3.56) 



This expression highlights both the non-perturbative nature of the expansion in terms 
of fractional powers of X/H as a result of the infrared divergences and the limit of validity 
of the leading order approximation determined by the pole in A for subhorizon modes: the 
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second term in the bracket begins to dominate for 



H 



> 



271 H 
A 



(3.57) 



Corollary and rules to extract the poles in A: From the explicit calculation in 
the appendix and the results obtained above we learn that the leading infrared behavior 
manifest as the leading poles in A correspond to the process of absorption and emission of 
superhorizon quanta depicted in fig. ([6]) for the self-energy for the cubic interaction. The 
pole arises from the integration in a small band of superhorizon wavevectors of width /i — )■ 
and can be extracted by implementing a set of simple rules that can be inferred from the 
example of the one-loop diagram corresponding to the self-energy in the theory with cubic 
interaction depicted in fig. (jl]) and given by eqn. f l3.34p . The momentum integral features 
two regions in which the infrared behavior leads to poles in A, p < yU and \k — p\ < fi, 
however, the second region is equivalent to the first by rerouting the external momentum, 
this yields an overall factor 2, namely the number of possibilities to reroute the external 
momenta along one internal lines. In the region of momentum integration p < fi we can set 
\k — p\ ~ k and extract the product gu{\k — p\,ri) g*{\k — p\,ri') ~ gu{k, rf) gl{k, i]') outside 
the integral which becomes 



zvr^ Jq Sit vr-^ yrji]' 



rji]' 



2A 



Sir'^Ariri' 



A 



\n[fi'^r]ri'] + 



(3.58) 



and the ln[/i] is cancelled by the integration with p > fi as explicitly shown in the appendix. 
Thus to leading order in A the self-energy (I3.34p becomes 



^io{k,ri,ri') 



^4^9:{k;v)9.{k;v')i2) 



H^ rji] 
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9u{k,r])gt{k,v') 
g^{k;v)\'\g,ik;i)\' 



Sir'^A'qri' 
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(3.59) 



the first term gl{k; ri)gt,{k; rj') corresponds to the external lines, the factor (2) is the number of 
lines through which the external momentum k can be rerouted, the first bracket corresponds 
to the internal line that carries the external momentum g^{k;r]) g*{k;r]'), the last bracket 
is the leading order contribution from the integration of superhorizon modes in the band 
< p < yU — )■ and the dots stand for subleading terms that do not feature poles in A. The 
expression (I3.59P confirms the leading order result of the self-energy (I3.44p . 

This result generalizes to the following set of rules valid for an irreducible self-energy 
graph with n + 1 internal lines and n loop integrations 

• Reroute the external momenta to run along one of the internal lines in the graph: there 
are n + 1 possibilities that give the overall factor n + 1, the "propagator" associated 
with this internal line is 

g,{k,v)g:{k,v'), (3.60) 



23 



For each loop momentum integral (with < p < ^) associate a factor 

1 



87r2 A (r/r/') 



(3.61) 



arising from the integration within a band of width /x — )■ of superhorizon quanta. 
Therefore the order of the pole in A is determined by the number of soft internal lines, 
hence the number of superhorizon quanta emitted by the decaying particle as depicted 
in fig. (p. 



The external lines correspond to 



9tik;v)9uik;r]') 



(3.62) 
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IpI < /".... ... \p\< f^ 



k k-p k k k-p k-p k 

(a) (b) 

FIG. 6: Processes that contribute to the leading order poles in A: (a) intermediate state of 
superhorizon modes, (b) emission and absorption of superhorizon quanta. 



D. Quartic self-interaction: 

We are now in position to use the above rules to obtain the leading order contribution to 
the self-consistent mass and decay width for the case of Ax^, whose interaction Hamiltonian 
is given by 



Hr 



d^x 



(3.63) 



2 if 2 ^2 
The one loop tadpole diagram shown in fig. (J3]) contributes to the diagonal matrix element 
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We obtain the leading order infrared pole of the two loop contribution by implementing 
the rules described above and summarized by the diagram shown in fig. ([7]), we find the 
leading order two-loops self energy to be 



^io{k,ri,ri') = 96 X'^ gl{k;r])g^{k; 7]') (3) g^{k,r]) gl{k,ri') 
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(3.65) 
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The dots stand for higher order terms in A that may contain logarithmic contributions, 
furthermore, just as in the case of cubic vertex there is a short distance divergence leading 
to both mass and wave function renormalization, these would have to be studied in detail 
and their assessment is beyond the scope of this article, however these are suppressed by a 
factor A^ oc A because the short distance divergences are the same as in Minkowski space 
time and of order A^ since these are impervious to the infrared divergences that lead to the 
poles in A. 



FIG. 7: Leading order in A for the two loop self-energy in A^^- The dashed lines correspond to 
the integration of the superhorizon modes within a band of width /x — )• 0. 



Again to leading order in A the two loops diagram yields a real contribution to S therefore 



no contribution to the mass. The double pole in A confirms the "infrared counting" of ref . |26 
but disagrees with the result of the "plane wave ansatz" in ref. [311] as in the case of the cubic 
coupling. Therefore the self-consistency condition becomes SM"^ = leading to the same 
result as the lowest order tadpole (12. 55112. 57p . 

Superhorizon modes: For —kr] ^ —krjo ^ 1 the rj- integrals yield the same result as 
for the cubic vertex using the result f l2.57p and to leading order we obtain in this case 



2A 



Ci,(r/) ~ e-^(-^'') ; 7(-A;r/) = ^ ^^ . (3.66) 



H 

kphiv) 



Note the power of A: the double pole in A leads to the result that the two loops contribution 
is of C(A) and the terms that have been neglected are at least of 0(A'^/^). As depicted in 
fig- (El) the decay process is dominated by the emission and absorption of two superhorizon 
quanta. Thus the power of A in the denominator is directly determined by the number of 
emitted superhorizon quanta. 

Furthermore, in analogy with the interpretation of the width of superhorizon modes for 
cubic coupling in its relation to the bi-spectrum, it is clear from fig, (ip that the imaginary 
part of S(A;, rj, r]') is related to the trispectrum of scalar fiuct nations |45|. and the power 1/k^ 
originates now on a highly squeezed configuration of the trispectrum]^ in which two of the 
momenta |p|, |^ < /i — )■ whereas the other two are ~ k. The integral over both bands 
of wavevectors yields the double pole in A, the order of the pole describes the number of 
superhorizon quanta emitted and absorbed and in this case the number of "squeezed" sides 
of the trispectrum configuration. 

Subhorizon modes: For —kr]Q ^ —kr] ^ 1 the calculation follows that in the previous 
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case, we find 
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(3.67) 



As in the case of tlie cubic coupling we expect tliat terms of (9(A^/^) multiply powers of 
kphiv)/H which will limit the validity of the leading order term in f l3.67p for modes deep 
inside the Hubble radius. We expect, just as in the previous case that for these deep 
subhorizon modes, the leading contribution will be determined by the conformal coupling 
limit, which clearly must be studied in detail for a deeper assessment, a task well beyond 
the scope of this article, which focuses on the long- wavelength limit. 



IV. CONCLUSIONS, COMMENTS AND FURTHER QUESTIONS 



A. Conclusions: 



Motivated by questions on the infrared stability of particle states in de Sitter space time, 
in this article we study massless minimally coupled scalar theories with cubic and quartic 
interactions focusing on infrared aspects and the time evolution of single particle states. 
In agreement with previous work [261. [ 



3l| we find that infrared divergences of massless 
minimally coupled fields lead to the self-consistent generation of a mass. In the case of 
the cubic coupling mass generation is a consequence of a radiatively induced expectation 
value of the scalar field which leads to the formation of a non-perturbative condensate 
{(f){x,ri)) oc H(^H/Xj as a result of the infrared divergences. This expectation value 



induces a self-consistent mass given by 



MocX'/'H^/' 



leading to a de Sitter invariant result 



')ocH' 



For a quartic coupling we find 



M oc A^/^ H 






^)oc 



Ml 



(4.1) 



(4.2) 



(4.3) 



respectively. The self-consistent mass generation results in that the infrared divergences 
in self-energies are now manifest as poles in A = M'^/3H^ as a result of the emission 
and absorption of superhorizon quanta. The two point correlation function of superhorizon 
fluctuations acquire an anomalous dimension 2A, and (0^(x, r^)) is de Sitter invariant. 

The self-consistent treatment of mass generation is combined with the non-perturbative 
Wigner-Weisskopf method introduced in ref. 40| to extract the time evolution of single par- 
ticle states. The lack of a global time-like Killing vector entails the lack of kinematic thresh- 
olds and that all single particle states decay into quanta of the same field. The decay is 
dominated by the emission and absorption of superhorizon quanta and hastens when the 
wave- vector of the single particle state crosses the Hubble radius. The radiatively generated 
mass regularizes the infrared behavior of the decay rate and leads to a rearrangement of 
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the perturbative expansion non-analytic in the couplings. We obtain a set of simple rules 
that yield the leading order contributions to the decay law of single particle states and we 
find that the amplitude of single particle states with superhorizon wavelength decays as 
Q-iM ■ ^(^^^ oc {X/HY^^ [H / kph{'r]yf' ; A [H / kph{riyf for cubic and quartic couplings respec- 
tively to leading order in M/H. The generation of mass as a consequence of the infrared 
divergences in radiative corrections results in a re-arrangement of the perturbative expansion 
non-analytic in the couplings. 

The power law 1/k^ in the decay "width" of superhorizon modes is expected on the basis 
of the relation between the self-energy and the bi-spectrum or tri-spectrum]^^ of scalar 
perturbations in the highly squeezed limit for cubic and quartic coupling respectively. The 
order of the pole in A describes the number of superhorizon quanta emitted and absorbed 
in the intermediate state and also the number of highly squeezed sides in the bi-spectrum or 
tri-spectrum configurations respectively. 

Mass generation and decay entails that all single particle states become quasiparticles 
during de Sitter inflation. 



B. Agreements and disagreements: 

In the case of quartic couphng Xcj)^, our self-consistent one-loop results (I2.56|2.58l) agree 



with those of refs.[26|, l3J] and the one- loop result i n[32 | and qualitatively with the same 
power of A with the general results obtained in refs. [25|, |29|, |32, |4J] but not quantitatively, 
the disagreement is in numerical factors of order one. We emphasize that the result for the 
radiatively generated mass fl2.56p does not receive contributions at two loops to leading order 
in A as explained in detail in section (jIIIDI) . The stochastic approach of refs. [25i. ,29] obtains 
a Fokker-Planck equation for a coarse grained average of the field and extracts (0^) from the 
asymptotically long time solution of this Fokker-Planck equation, in ref. [29] such solution 
features a factor r(3/4)/r(l/4) with respect to the one-loop result, a different numerical 
factor is found in ref.|25|. It is not clear from the stochastic approach which type of diagrams 
are being resummed by the Fokker-Planck equation, in particular, whereas the two-loop 
diagram studied in section (lIIIPp is non-local (and retarded), the stochastic approach does 
not seem to reflect any retardation or non-locality effects. The coarse grained field is treated 
as a classical stochastic variable whereas in a diagrammatic approach the loop integrals 
reflect the contribution of intermediate states and it is by no means clear (at least to this 



author) how such contributions contribute to the Fokker-Planck equation (although ref. [30| 
provides a field theoretical justification, precisely how the higher loop retarded contributions 
are summed up is, again, not clear to this author). These numerical factors are also obtained 
in the result of ref. [441] which presents a completely different approach: here a Euclidean 
formulation with compact time on the sphere is implemented, the zero mode (on the sphere) 
is isolated and the path integral of one single mode is carried out yielding the result (0^) oc 
if^/vA with numerical factors similar to those of the stochastic approach. However, it is 
by no means clear (at least to this author) how the Euclidean formulation captures any 
dynamical information which is manifest in the stochastic approach, perhaps the Euclidean 
compactified description captures the asymptotically long time solution. Furthermore, in 
the Euclidean compactified formulation, the zero mode is separated from the other modes 
by a gap and can be isolated unambiguously, whereas in the Lorentzian version the poles 
in A arise from the integration of a band of infrared wavevectors. Finally, whereas the non- 



perturbative results beyond one loop obtained in|32| purport to provide a full resummation 
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of loop diagrams, the final result (eqn. (55) in|32[) also disagrees quantitatively with the 
results from the stochastic approach [25|, l29| and the compactified Euclidean |44i] approach, 
furthermore, the plane wave ansatze proposed in ref. |32| is manifestly a local approximation, 
and the authors state clearly that they assume that the infrared effects can be resummed in 
a (local) mass term. This is not borne out in the results obtained in section (IIIIPp which 
reveal that infrared effects also affect the decay of quantum states and not only the mass. 

An important aspect that transpires from the two loop result in section (lIIIPp is that 
the limit fc — )■ is singular, and in fact in separating the integration band g < /i — > we 
have manifestly kept the external momentum k ^ 0. It is conceivable that the solution to 
the discrepancy lies in a more thorough treatment of the k = limit, after all both the 
stochastic approach and the Euclidean approach consider only the "zero mode". But if 
this is the case there seems to be a discontinuity in the treatment of this mode. Perhaps 
treating this single mode as a condensate (but without manifestly breakirig the underlying 



symmetry) and separating its contribution as seemingly advocated in ref.[l8j is the correct 
procedure. 

In summary there is general qualitative agreement among the various different approaches 
that in the case of Xcp'^ there is a radiatively generated mass M oc X^/^H arising from the non- 
perturbative build-up of infrared effects. However there is a quantitative disagreement in the 
proportionality constant by numerical factors of 0{1). The origin of the discrepancy is diffi- 
cult to extract because the different approaches that purport to provide a non-perturbative 
resummation cannot be directly compared, while the stochastic and Euclidean approach 
share similar numerical factors in the corresponding results, a comparison between the two 
approaches and with the diagrammatic approach is far from obvious or clear. Understanding 
the origin of this quantitative discrepancy is certainly worthy of study, but clearly beyond 
the purview of this work. 



C. Comments and further questions: 

• Powerful results on the asymptotic behavior of correlation functions in Euclidean de 



Sitter were obtained in references [3 71. |4J, |47||. In particular for heavy fields with bare 
mass M^ > 9H'^/4 the decay of correlation functions was interpreted as an imaginary 



mass, a result confirmed in ref.[38|. The continuation from heavy mass to the light 
mass (or massless) case considered in this article is neither direct nor clear as the 
superhorizon limit seems to be different, as exhibited in the integral in eqn. (I3.5ip 



which is also found in ref. [111 . I38l | . The analytic continuation from Euclidean results 
to the Lorentzian de Sitter case if valid also for the massless or light mass case would 
yield a powerful method to extract the decay of correlators and perhaps establish an 
equivalence with the decay of quantum states. These aspects clearly merit further and 
deeper study. 

A non-perturbative resummation of infrared divergences should yield correlation func- 
tions that are well behaved in the superhorizon limit and as ri — )■ 0~. While resum- 



mation methods such as the dynamical renormalization group [261. I27j may ultimately 
be a successful approach, in this article we followed a different route, namely to study 
directly the time evolution of states adapting and generalizing a resummation method 
that has proven successful in other areas of non-equilibrium phenomena. In ref. [40| the 
equivalence of this method with the dynamical renormalization group was established 
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in Minkowski space time and it would be fruitful to establish a similar relation in de 
Sitter cosmology. The main strategy that we advocate can be best described by the 
example of the two point correlation function of a pion field in Minkowski space time. 
The pion decays into a lepton pair, the self-energy features a two lepton threshold 
and the retarded propagator features a complex pole in the second Riemann sheet as 
befits a resonance. The time evolution obtained from the frequency Fourier transform 
reveals the exponential decay of the correlation function oc e~'"*/^ with T the decay 
width. Now consider the pion correlator (7r(fc, t) 7r(— fc, 0)) introducing a complete set 
of states between the two fields. Passing to the interaction picture, the one-pion in- 
termediate state is the solution of the Wigner-Weisskopf equation f l3.15|3.16p with a 
two lepton intermediate state[425, and decays as e~'"*/^ , which is the correct long-time 
limit from the time Fourier transform of the full propagator. The Wigner-Weisskopf 
solution is equivalent to a Breit-Wigner approximation to the full propagator [40, 42|. 
We will report on this approach to obtain the correlation functions as generalized to 
de Sitter space time in a future study. 

The decay of single particle excitations into superhorizon quanta hastens as the 
wavevector crosses the Hubble radius. This process entails the creation of particles 
leading to a build-up of the population of the produced particles which may affect the 
time evolution of the single particle states as the produced particles recombine into 
the initial state. This possibility requires to study a Boltzmann equation in which 
both the decay process 1 — )■ 2 and its reverse, the recombination 2 — )■ 1 are taken 
into account. The question to study is whether a detailed balance emerges where the 
recombination process balances the decay reaching a steady (or perhaps equilibrium) 
state. An assessment of these processes requires to obtain a Boltzmann equation as 



reported in refs. [48|, |49[ for the case of heavy fields, but adapted to the massless case 
and including the self-consistent mass generation mechanism as in the case studied in 
ref. J50| in Minkowski space-time. This program is relegated to future study. 

• The relation between the decay "width" of superhorizon modes and the bispectrum or 
trispectrum of scalar fluctuations raises an interesting question: are non-gaussianities 
of curvature perturbations |45| related in any way to the decay of either adiabatic 
or isocurvature superhorizon fluctuations?, if so is there any imprint on the cosmic 
microwave background anisotropics?. Clearly this question also merits further study. 
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Appendix A: Calculation of the Self Energy (|3.38|) 
Consider the integral in f l3.38p 

/■oo 

I[k,ri,ri']= 'Hu{p;ri',v') J^u{k + p; vv') ~ J^u{\k - P\i VV) P dp (Al) 

Jo L -I 
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introduce an infrared cutoff /i — )■ and write 

I4k,r],r]'] = Io + Ii-l2 
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Obviously the total integral I[k, rj, rj'] is independent of the cutoff /i. Since for p — > oo 

'H,{p-r]-r]')(xe-'P^'>-'>''^ (A6) 

we introduce a convergence factor 

rj — rj' ^ rj — rj' — ie ; e — )■ 0"*" . (A7) 

The integral Ji is infrared finite as long as /i 7^ and to leading order in A we can set 
z/ = 3/2, namely A = therefore 
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In the integral I2 we must isolate the region p c:^ k extracting the pole in A and outside this 
region we can replace z/ = 3/2; A = 0. Therefore we write 
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In 1 2 and 1 2 we can set z/ = 3/2; A = since these are infrared finite, namely 
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For I2 and I2 we find 
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therefore keeping the pole in A and the leading logarithmic infrared contribution as /i — )■ 
we find 



rib) Ac) 
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rjrj' 



where 



Furthermore we find 
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The integrals are straightforward, we find to leading order in A 
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where Ci\ Si are the cosine and sine integral functions respectively and 
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with 7 the Euler-Mascheroni constant. The infrared logarithm ln[/i] cancels out as it can be 
easily seen from the expressions (JASHASP by taking derivatives with respect to /i. 

To write the first term in (]A17p we have used the definition of T>[k, rj, rj'] eq. (1A14P and 
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The second term in the first line in (IA17P is recognized as the contribution from the 
massless conformally coupled case|28|, |40|, and can be written as 
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where e is a physical short distance cutoff that reffects the logarithmic ultraviolet divergence. 
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